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Abstract 

We analyze Nf — 2 nucleon mass data with respect to their dependence on the pion 
mass down to = 157 Me V and compare it with predictions from covariant baryon 
chiral perturbation theory (BChPT). A novel feature of our approach is that we fit 
the nucleon mass data simultaneously with the directly obtained pion-nucleon cr-term. 
Our lattice data below = 435 MeV is well described by 0(p 4 ) BChPT and we find 
a = 37(8) (6) MeV for the cr-term at the physical point. Using the nucleon mass to set 
the scale we obtain a Sommer parameter of r = 0.501(10)(11) fm. 

Keywords: nucleon mass, pion-nucleon sigma term, Sommer scale, covariant baryon 
chiral perturbation theory, finite size corrections 



1. Introduction 

Predicting low-energy hadronic properties is one basic goal of lattice QCD. A partic- 
ular challenge to experiment as well as to theory is posed by the so-called pion-nucleon 
cr-term 

a = m e (N\(uu + dd)\N) (1) 

which parametrizes the light quark contribution to the nucleon mass. Here mi = m u = 
■nid denotes the light quark mass. 
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At present, phenomenology does not give a clear picture of the magnitude of a. A dis- 
persion theoretical analysis led to a = 64(8) MeV [I], similar to the value a = 64(7) MeV 
obtained later in Ref. [2J. However the analysis of Ref. [3] suggested a much lower value, 
a = 45(8) MeV, which was also found in Recently, a new evaluation resulted in 
a = 59(7) MeV [5]. 

Calculating the pion-nucleon cr-term directly on the lattice is a computationally inten- 
sive task as it involves the computation of quark-line disconnected correlation functions. 
This has become feasible recently and was performed by us and others [SHE] > though only 
at a single value of ni£ and of the lattice spacing a. An alternative, which has often been 
used in the past (see, e.g., the recent studies [8Ullj). is given by the Feynman-Hcllmann 
theorem. It allows us to express a in terms of the derivative of the nucleon mass Mjv 
with respect to mf. 

9M N 

Calculating the nucleon mass on the lattice as a function of mi is relatively straight- 
forward, but becomes expensive close to the physical point. 

Until very recently, lattice QCD calculations have therefore been performed at rather 
large quark masses, such that results had to be extrapolated to the physical point over a 
wide range. Thanks to the efforts of different lattice QCD collaborations during the last 
years this situation has much improved. The QCDSF Collaboration, for example, has 
generated a large set of Nj = 2 gauge field configurations for a variety of quark masses, 
lattice spacings and volumes, reaching down to pion masses of about 157 MeV. These 
simulations employ the standard Wilson gauge action and the non-perturbatively O(a) 
improved clover action for two flavors of mass-degenerate quarks. 

In this paper we analyze nucleon mass data obtained from these simulations with 
respect to their quark-mass and volume dependence, compare this to SU (2) covariant 
baryon chiral perturbation theory (BChPT), and extract a value for the pion-nucleon 
cr-term utilizing the Feynman-Hellmann theorem. A novel feature of our analysis is 
that we combine the nucleon mass data with a direct determination of the pion-nucleon 
cr-term [7], fitting both simultaneously to the corresponding 0(p 4 ) BChPT expressions. 
It turns out that this gives much more reliable and precise results compared to fits to 
the nucleon mass data only. 

The outline of this article is as follows: In the next section we summarize and describe 
our lattice data. Our fitting procedure is described in Sec.|3l The results of the fits 
are then discussed in Sec. [4] In Sec. [5] we summarize the outcome of our analysis. In 
|Appendix A| we review the results from BChPT which underlie our analysis. Some more 
details on the fits can be found in Appendix B 



2. Lattice data 

Our lattice data for the pseudoscalar (am^) and nucleon mass (aAfjv) were extracted 
from one-exponential fits to smeared-smeared correlators. For the smearing we used 
Jacobi smearing for all data sets apart from the more recent analyses at = 5.29, 
K = 0.13632 and k — 0.13640. For k — 0.13632 our aforementioned direct calculation 
of the cr-term was performed [TJ, for which a more optimized smearing was critical to 
obtaining a signal for the scalar matrix element. We achieved this by using Wuppertal 
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Table 1: Lattice data for the pseudoscalar (am n ) and the nucleon mass (aMjv) in lattice units. In 
columns 6 and 7 we list the corresponding values in units of ro. For am w and aMjv we give the statistical 
errors, for r^m^ (ro-Mjv) the errors are the combined statistical errors of rrj/a and am n (aMjv). Stars 
in the first column mark all those entries which enter our fits [depending on the upper limit set for 
(r m„) 2 ]. 



smearing with APE smoothed links. This improved method was also applied to our 
calculations at k = 0.13640. 

TableQ] lists our data for the pseudoscalar (am^) and nucleon masses (aMjv). We 
also state the corresponding values r^m^ and roMjv in units of the Sommer scale ro, 
extrapolated to mi — [T^] (see Tablc[2|. The error quoted for am^ and oMjv is the 
statistical uncertainty of the data, while for r^m^ (ro-M/v) it is the combined error of 
ro I a and am^ (clAIn). 

In Fig. [l]we show our data for roMjv plotted versus (rom T ) 2 . Different symbols and 
colors are used to distinguish between data for the different f3. Data points which refer to 
the same (/3, n) but different lattice volumes are connected by dotted lines to emphasize 
finite volume effects. 

From this figure (and also from closer inspection of the data) we find that, within the 
given precision, our data show no systematic dependence on the lattice spacing and, thus, 
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Figure 1: Lattice data for roMpf versus (rom^) 2 for different lattice spacings (/? = 5.25, 5.29, 5.40) and 
volumes. Open (filled) symbols refer to points where L < 3ro (L > 3ro). A black-framed (yellow) circle 
indicates the physical point assuming ro = 0.5 fm. The dotted lines connect points of the same (0,k) 
but different lattice size. They are meant to guide the eye and to illustrate finite-volume effects. 



are close to the continuum limit. We therefore scale our data by the respective values 
of the chirally extrapolated ro/a and do not attempt a continuum- limit extrapolation. 
Finite volume effects, on the other hand, are clearly visible and will be incorporated 
into the fits. In fact, they provide an additional valuable input, because some of the 
low-energy constants (LECs) also enter the volume corrections. 

Before fitting these data to BChPT expressions, all values for the pion mass have to be 
extrapolated to infinite volume. A method to calculate these finite- volume corrections has 
been worked out in |13j . When applying this method to our data we find, however, that 
the corresponding values for m n L have to satisfy at least m^L > 3.5. Below this limit, 
the finite-size effects of our pion mass data are stronger than the calculated corrections. 
If m n L > 4 the data points and the extrapolated values even agree within errors (see 
Fig. [2] for an illustration). This means that, within the available precision, correcting for 
finite-volume effects according to [TS] will only help us if 3.5 < m n L < 4.0. In all these 
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Table 2: Lattice estimates of ro/a for different /3, extrapolated to the chiral limit 1121 . 
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Figure 2: Volume dependence of our pion mass data. Left top: /3 = 5.29, k = 0.13632. Left bottom: 
P = 5.29, k = 0.13640. Right top: = 5.40, re = 0.13640. Right bottom: = 5.40, n = 0.13660. 
Dashed-dotted lines represent finite-volume extrapolations, obtained from applying the method of |13| 
to the point for the smallest ro/L. Vertical dotted lines denote constant m^L. Gray areas mark the 
region (m^L < 3.5) where finite-volume corrections are not under control. 



cases we have at least one point with m^L ps 4 (or larger). 

For the final fits to BChPT, we therefore exclude all (/?, k) combinations for which 
there is not at least one data point that satisfies m^L > 3.5, and then take the value for 
r^m^ from the largest available volume (with m n L > 4) as our estimate of the infinite- 
volume limit. As argued above, a finite-volume correction according to Ref. [12] would 
not give different results for rom T , but in this way we are a bit more conservative about 
the error. 

There is, however, one data point where we allow for an exception to this rule: the 
estimate for r^m^ at j3 = 5.29 and k = 0.13640 (m^ w 157MeV). At these parameters 
the largest available lattice size is currently 48 3 x 64, but we expect no drastic changes of 
rom % if the volume is enlarged further, because the PCAC-mass dependence of our pion 
mass data below 290 Me V (considering largest-volume data only) is quite linear already]^] 

As mentioned above, a novel feature of our analysis is that we take additional data 
for a into account, which comes from a direct determination. So far, we have computed 
this at j3 = 5.29, k = 0.13632 on a 32 3 x 64 and 40 3 x 64 lattice [7]. This corresponds to 
a pion mass of about 290 MeV, and we will use the value 

r a = 0.273(25) (3) 

from the 40 3 x 64 lattice in what follows. 



1 Note also that if with larger volumes the pion mass further decreased the Xr _vame °f our tits presented 
below would actually improve, with negligible effects on the fitted parameters. We have tested that. 
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3. Fitting to covariant chiral perturbation theory 



3.1. Fitting formulae 

We intend to fit our data t o BChPT. The expressions at next-to-leading one-loop 
order are given in Appendix A We expand those up to order ra 4 .: 
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These expressions involve the low-energy constants c\ , c 2 and C3 as well as the renormal- 
ized counterterm coefficient e\ and 
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which depend on the renormalization scale A. The pion decay constant and the 
nucleon axial coupling constant gA are taken at the physical point, which is consistent 
with the order of BChPT we are using. For the fits we will set F^ — 92.4 MeV and 
g A = 1-256. The renormalization scale A is set to A = mP hys = 138 MeV. f\ 

To correct for finite-volume effects in our nucleon mass data we employ the finite- 
volume correction at next-to-leading one-loop order in BChPT. It reads 

AM N (ml,L) = AM^(ml,L) + AAf( 4 )(m 2 ,L) , (7) 

where AM^ and AM' 4 ' can be read off from Eqs. (IA.15I) and (|A.16| with the substi- 



tution m — > m^. This is consistent with the order we are using. 

Since our data is given in units of r , for the fits we have to re-express all dimensionful 
quantities in units of ro, too: 

Mjv^tt) — > Mjv(TO,r), AMjv(m 7r , L) — > AM/v (to^ , L) and o^m^) — > (r(m w ) (8) 

where a indicates the expression is understood in units of ro, e.g., Mjq = tqMn. 

As not all parameters of Eqs. Q and ^ are well constrained by our fits, some of 
these will be fixed to their phenomenological values, e.g., F„, c 2 , C3 and ^3. Consequently, 
a physical input value for ro has to be provided as well, which is however unknown a 
priori. 



2 Note that choosing another renormalization scale will only change the values of and ij. All other 
parameters are invariant. We have checked that this is satisfied for our fits. 
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We fix r for each fit separately by iterating over different physical (k = 1, 2, . . .), 
plugged into the fitting formulae until ro comes out self-consistently from the fit, that is 

e > r< fc) - r * X) where r« = ^— ^ >- . (9) 

M N 

We set e = 0.001, and for most of our fits ro converges after 10 to 15 iterations. By 
construction, all fits will pass through the physical point Alf^ ys = 938 MeV at mP hys = 
138 MeV. 

Two kinds of fits will be discussed below: fits to our nucleon mass data, using the 
% 2 -function 

N d 



W (M N (xi) + AM N (xi, yi ) - ZiJ 



xl=Y. ^ 72 (io) 

and combined (simultaneous) fits to the nucleon and cr-term data, using the x 2 -function 

xL = ti+ Mxi l-~ Zi)3 ■ (ID 

Here x = r^m^, y = L/r , z = r M^ refer to our measured points and e denotes the 
error of r^M^. In Xat<t' %j (j € {1, . . . , Adata}) refers to our single directly determined 
result for ro<7 and e-j is its error, see Eq. As finite-size effects appear to be negligible 
for this number (see Ref. [7]) we do not apply any finite- volume corrections in this case. 

3.2. Fit ranges 

Our nucleon mass data of Table [l] covers a range of rom v values from 0.42 up 
to 4.04, and L/tq ranges from 1.71 to 6.85. In physical units this corresponds to 
m w w 0.17... 1.58 GeV and L = 0.85... 3.4 fm when r = 0.5 fm. When fitting these 
data we do not know a priori for which range of values of and L we can trust our 
fitting functions. We therefore vary constraints on r^m^ (£/r*o) from above (below). 
The constraint L/vq > 3 turns out to be low enough such that there is sufficient data to 
perform stable fits but also large enough so that AMn captures the finite- volume effects, 
see, e.g., the discussion below. 

For r m w , on the other hand, we find it reasonable to constrain it from above by 
{Tom^y^^y. = 1.6, if not (rom 7r ) 2 lax = 1.3. Both these upper bounds give results which 
agree within errors, albeit with a larger uncertainty for the latter. Also our independent 
measurement of r$o [Eq. ([3])] at TQm„ w 0.735 is then well reproduced by Eq. ([5]), not 
only for a combined fit [Eq. ( [TT] )] but also if one uses Eq. (J5j) with the parameters from 
a stand-alone fit to the nucleon mass data. For larger (rom^j^, say (r m T )^ lax = 3.0, 
the fits change quantitatively and qualitatively: The m^-dependence of the nucleon mass 
[Eq. Q] becomes more concave in shape and our data point for r^a lies below the fit 
curves. If one only considered the nucleon mass data, such fits over a larger m^-range 
would roughly capture the overall m^-dependence [even up to ss 1 GeV where 0(p 4 ) 
BChPT certainly does not hold], but this is likely to be accidental for the given order 
as was pointed out already in |14j . Our definitive analysis will therefore be restricted to 
fits for which (ro ?7i 7r ) max ^ 1.6. For our final estimates of tq and o we will even restrict 
ourselves to (roTOTi-)max < 1.3. 

7 



3.3. Parameters 



Let us now discuss the parameters in the fit functions. Besides F v and gA, the 



functions in Eqs. (jlOh and (11) contain five and six free parameters, respectively. These 



are M , Q {i = 1,2,3) and e\ for Eq. (10) and Z3 in addition for Eq. (11). Ideally one 
would like to determine these all from fits to the lattice data. This is however not possible 
with our current data and we therefore have to fix some of the parameters to values from 
the literature. 

As I3 only enters r$o [see Eq. ^} this parameter should not be left free in a fit with 



only one data point to constrain it. For our combined fits [Eq. ( 11 )] we therefore fix it to 
the FLAG-estimate [TS] 

l 3 = 3.2(8) (12) 

and check for the stability of our final results varying I3 = 3.2 within one standard 
deviation. 

For C2 and C3 we proceed similarly. At first one might be tempted not to fix these two 
parameters at all. However, like the renormalization-scale dependent parameter e£, C2 
and C3 gain influence at larger m w . So, if these parameters are all left free in fits to low 
m„-mass data, their uncertainties would be unreasonably large. We therefore decided to 
fix C2 and C3 to their latest phenomenological values [TH [17] , 



c 2 = 3.3(2) GeV -1 and c 3 = -4.7(1.3) GeV" 1 , (13) 

and investigate the stability of our final results varying the less precisely known parameter 
C3 by one standard deviation. Additional fits are performed where C3 is left as a free 
parameter. 

To check the stability of the fit parameters c\ and Mq, and also of our estimates of 
roAfJ hys and Tq<j at the physical point, additional fits are performed where the expressions 
in Eqs.Q and ^ are truncated at orders O(m^) and 0{m%). 

4. Results 

4-1. Discussion of our fits 

In Fig. [3] we show an example fit to the nucleon mass data. Full (black) circles in 
Fig. [3] represent our lattice data for the nucleon mass, and the surface is a fit to these 
points. We fit to the surface ro-Mjv^o^Vj L/tq), hence the finite-volume corrections are 
determined directly through the fit as well. The (red) diamonds in Fig. [3] represent the 
lattice data, after subtracting these volume corrections. 

An overview of our fits and parameters can be found in | Appendix B| where we list all 
our combined and stand-alone fits in Tables |B.4| and |B.5| respectively, and also provide 
more detail for the interested reader. 

In Figs. |4] and [5] we display the fit curves corresponding to some of the results listed 
in these tables, together with the lattice data after subtracting the volume corrections. 
The overlap of points indicates the quality of these fitted corrections. In Fig.|4j three 
of our combined fits are shown (Soo3, Soo2 and Sool of Table lB~4|) . each for the same 



choice of fixed parameters (02, C3 and I3 as given in Eqs. (12) and (13)), but for different 
fit ranges 

(t-oto^) 2 < (rom^)^ , (14) 




Figure 3: Fit (surface) to the nucleon mass data (full circles) for a range of ro/L and (rom T ) 2 values. 
Black circles lie above the surface, gray (or half gray) below (or on) the surface. The line and the red 
diamonds at ro/L = mark the fitted infinite volume prediction. For the sake of simplicity, for each 
(fO m ir) 2 omv one extrapolated (red) point is shown at ro/L = 0. 



where ('"o m ir)max * s 3-0, 1.6 or 1.3, always requiring L/r a > 3. Fig.[5]shows a correspond- 
ing stand-alone fit (Noo2 of Table lB~5t to the nucleon mass data for (j"o"i 7r )^ lax = 1.6. 

Our fits perform significantly better if the er-term constraint [Eq. ([3|] is incorpo- 
rated. For these combined fits the individual uncertainty of each fit parameter is smaller 
than for the corresponding fit to the nucleon mass data alone, while the Xr~ va l ues 
(Xr = X 2 / n( ¥) largely remain unaffected. Fits qualities are inferior, if nucleon mass 
data up to (rom n ) 2 = 3.0 is included, but if one lowers this upper bound to 1.6 or 1.3, 
X^-values around 1 can be reached. Also the data point for r$o [Eq. ^] then agrees with 
the 0(p 4 ) BChPT expression for cr(m T ) [Eq. JH])], whether or not the data point for r^a 
was included in the fit (compare Figs. H and pi. For the fit range (rom^) 2 < 3.0 this is 
not the case anymore. 



4-2. Weighted averages 

Based on the results summarized in Table [B~4l we estimate the weighted averages of 
our fit results for ro<7 p hys and roMf^ ys . For the weights we use the statistical error, and 
only values from fits with xf- < 1-3 are allowed to enter the average. We obtain the values 
listed in Tabled As one can see from this table, our results for tq for (rom,,-) 2 < 1.6 and 
(rom-x) 2 < 1.3 are consistent within errors, including that from the stand-alone fit. 

In this table we also give the systematic error due to varying c 3 (second parenthesis) 
and I3 (third parenthesis) one standard deviation around their phenomenological values 
in Eqs. (Tl2|) and (13 1. In total the systematic error is as large as the statistical error. 
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Figure 4: Simultaneous fits to the nucleon mass (left) and tr-term data (right) for three fitting windows 
with (^0"i7r)max = 3.0, f.6 and f.3 (from top to bottom). These fits are labeled Soo3, Soo2 and Sool 
in Table^O] where C2 = 3.3GeV _1 , C3 = — 4.7GeV — 1 and Z3 = 3.2. Lines, error bands and (full red) 
points are shown for the limit L — ¥ 00 (cf. Fig.[3]|. A black- framed circle marks the location of the 
physical point using — for each plot separately — the ro-value for which the fit is self-consistent. Open 
points did not enter any fit. In the left plots, the overlap of red points at (rom^) 2 = 0.436 and 0.538 
indicates the quality of the (fitted) finite- volume corrections. 
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Figure 5: Left panel: stand-alone fit to the nucleon mass data (see fit Noo2 in Table lB~5t , Right panel: 
comparison of the cr-term data (open symbols) for a 32 3 X 64 and 40 3 X 64 lattice [J] and the BChPT 
expression for cr(m„-) with the parameters Mo, ci and e}. taken from the nucleon mass fit shown in 
the left panel. The black-framed circle is the value for tqM^ (fofjv) at the physical point for these 
parameters. In both panels, the one-sigma error band is shown in gray. As in Fig. [4] the overlap of (red) 
points (left panel) at same (rom^) 2 indicates the quality of the (fitted) finite-volume corrections. 



Fit 




(n)WW)max 


FoCphys 


r [GeV" 1 ] 


r [fm] 


a phys [MeV] 


oK) 


Xn 


1.6 


0.103(23)(5) 


2.51(6)(2) 


0.495(12) (4) 


41(9) (2) 


0(mt) 


2 

XNcr 


1.6 


0.095(11)(12)(4) 


2.54(3)(4)(2) 


0.501(6) (8) (4) 


37(4)(5)(2) 


0(mt) 


2 

XNcr 


1.3 


0.093(20)(15)(5) 


2.54(5)(5)(2) 


0.501 (10) (10) (4) 


37(8)(6)(2) 


0{ml) 


XjVct 


1.0 


0.121(5) 


2.49(3) 


0.491(6) 


49(2) 


0(ml) 


2 

XNcr 


1.0 


0.065(7) 


2.58(3) 


0.509(5) 


25(3) 



Table 3: Weighted averages of fit results with \ 2 < 1.3. The first row gives averages for (stand-alone) fits 
to the nucleon mass data; the remaining rows for (combined) fits to the nucleon mass and (T-term data. 
The first column specifies the order of the chiral expansion, the second the upper limit on r^m^. Some 
of the weighted averages come with a statistical and systematic error: The error in the first parenthesis 
is always the statistical error and the second (third) the systematic error estimated by changing the fixed 
parameter eg (T3) by one standard deviation [see Eq. ^12\ and j!3) ]. The numbers that went into the 
averages are listed in Tables [B.4| and |B.5| 



4-3. Fits to lower order expansions 

It is interesting to check the robustness of the above estimates for r Q and cr p hy S by 
fits to 0(m%) and O(m^) BChPT. Up to these orders, only c\ and Mq are left as free 
parameters. Note that we still have to correct for the finite-volume effect in the nucleon 
mass data. We do this by setting (as above) C2 = 3.3 GeV" 1 and C3 = — 4.7GeV _1 in 
AM N . 

For the fits to O(m^) and 0(771^.) BChPT we employ the x 2 -function for our combined 
fits [Eq. ( 11 )] . The fitting ranges are chosen as above but we add to these (f o^^Omax = 1-0- 
It turns out that only for (rom^)^^ — 1.0 reasonable fits to O(m^) and 0(m%) BChPT 
can be found. 

Our results for r-o<7 p hys and r$ from these fits are listed in Table[3] To ease the 
comparison we also show them with our 0(m\) results in Fig. [6} Open (full) symbols 
correspond to fits where > 2 (xl < 2), black- framed full symbols represent good fits 
where Xr < 1-3- 
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Figure 6: Results for rga at the physical point (left panel) and ro (right panel) from fits to BChPT 
truncating the fit functions M^f(m n ) and a(m n ) at different orders in m^. Only results for our combined 
fits [Eq. poll] are shown. Points are grouped together with respect to the upper cutoff on (rom n ) . 
Red diamonds are for fits to order O(m^), blue squares are for 0(m^). Green circles refer to the 
weighted averages in Table[3] where for the error the statistical and systematic errors have been added 
in quadrature. Full (open) symbols refer to fits for which \r < 2.0 (> 2.0); black-framed full symbols to 
fits where \r < 1-3. The gray bands represent the errors for the lower green circle, which corresponds 
to our final values for ro0" p h ys (left) and ro (right). 



As can been seen from this figure, fits to different orders in result in slightly 
different estimates both for ro and rocxphys, and also come with a different fit quality. 
However, these deviations get smaller when increasing the order or decreasing (romr)^. 
For fixed (^o^7r)max' points for different orders seem to alternate around a yet unknown 
value with the tendency of coming closer to that with each order. Results from fits to 
even larger orders or lower (fo TO 7r)max wm hkely be found within the error bounds of the 
0(m£) results. 

We therefore conclude that our estimates for ro and cr p hys, that is, the weighted 
averages of our results from the combined 0(m^) fits with (?"o?Ti 7r )^ lax = 1.3, lead to 
sufficiently conservative errors to accommodate all of the uncertainties involved when 
fitting nucleon mass data to BChPT. 

4-4- Comparing different orders 

Let us finally try to get an impression of the convergence properties of the BChPT 
formulae we are using. To this end we plot in Fig. [7] a combined fit (Sf ol in Tablc lB~4l 
where C3 is a free parameter) along with the curves which result from truncating the 
BChPT function at O(m^) and at O(m^), using the same parameter values. In addition 
we show two (three) curves where the contribution of fifth order in m w has been added to 
the fitted function varying the new LECs d\§, d\% and l\ appearing in this contribution 
within a phenomenologically acceptable (slightly expanded) range. At this order, the 



LECs dig and d\% enter as the difference 2d\ & — d r l& (see Appendix A for details) , which 
currently is known only approximate, 2e^ 6 — d\ & — (—2.0 ± 2.5)GeV -2 [T51 ITU] , For l\ 
we use the value for Z4 given in [15] . 

Given this range of expected values, we see that the fifth order correction becomes 
a non-negligible effect already at pion masses well below the physical kaon mass. More 

12 



138 200 300 400 m„[MeV] 
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Figure 7: Comparison of one of our combined fits to the contributions from lower orders (O(m^) and 
0(m^.)) and to the 0(m^.) expansion. For all expansions the parameters are fixed to those of the fit 
(Sf ol in Table lrTZll . The additional parameters for the 0(m|) expansion were fixed to values consistent 
with current expectations (see text). In addition, we show the O(m^) expansion assuming a larger 
value for 2die — dig than currently expected from phenomenology (very right dashed-dotted line). For 
simplicity, we show only finite-volume corrected data (diamonds) from the largest lattice volumes, even 
though for (rorrt,r) 2 = 0.435 and 0.540 also points from smaller volumes entered the fit. 

specifically, at pion masses of ~ 350 MeV, the fifth order contribution is already of about 
the same size as the third order term from the leading-one-loop correction. This leads 
us to the conclusion that, in case of the nucleon mass and sigma term, BChPT (with 
the current LECs) shows no sign of convergence beyond > 250 MeV. Discussions 
pointing in this direction can also be found in, e.g., [2UH2H] ■ If one allowed, however, 
for a slightly wider range for 2d\§ — c^ 8 , say 2d\ & — d\% — 3.0 GeV -2 , the fifth order 
contribution could be much smaller (see Fig. [7]). But this is speculative only, and an 
improved knowledge on the values for d[ 6 and d[ 8 is required to decide on that. At least 
from the small difference of the 0(p 3 ) and 0(p 4 ) functions at small (rom OT ) 2 , we learn 
that the LECs C2 , C3 and e\ are not well constrained by nucleon mass data at pion masses 
m„ < 300 MeV. 

5. Conclusions 

We have presented Nf = 2 QCD nucleon mass data. The corresponding pion mass 
values range from about 1.5 GeV down to 157 MeV. To estimate the nucleon er-term we 
have performed two kinds of fits to expressions from 0(p 4 ) BChPT: (stand-alone) fits to 
our nucleon mass data and simultaneous fits to the nucleon mass and cr-term data. The 
latter was determined in a separate study [7] at a pion mass of about 290 MeV. 
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Figure 8: Comparison of lattice estimates for the pion-nucleon cr-term at the physical point for Nt = 
and Nf = 2 + 1 QCD I9l lllll24"rl29| . The lowermost point represents our estimate for <r p hys [Eq. l|16[l] 



For the fits, different fitting ranges in m n and L (spatial lattice extension) have 
been tested. We find that if one demands m, < 500 MeV [(rom^) 2 < 1.6] and also 
L > 1.5 fm acceptable fits to 0(p 4 ) BChPT can be foundj^] These fits do not only give a 
good description of the dependence of the data but also of the finite-volume effects. 
Generally, our simultaneous fits perform better than fits to the nucleon mass data alone. 
They are also robust against variations of (? , o TO 7r) 2 n ax an d of the low energy constants 
C2, C3 and Z3. We have found a strong correlation between the counterterm coefficient e\ 
and C3. More data points below 500 MeV pion mass will be needed to resolve this issue 
or to fix C2 and C3 through lattice data. 

As our final estimates we quote 

r = 0.501(10) (11) fm (15) 

and 

cr phys = 37(8)(6)MeV. (16) 

These numbers are the weighted averages given in the last row of TableEl adding the two 
systematic errors in quadrature. These numbers result from our fits to 0(p 4 ) BChPT, 
fitting simultaneously the nucleon mass and cr-term data up to pion masses of about 
433 MeV. Within errors, we find these numbers to be consistent with corresponding 
fits to 0(p 2 ) and 0(p 3 ) BChPT (see Fig. [6]), and also with our recent estimate cr p hy S = 
31(3)(4)MeVfor N f =2 + 1 [9]. 



3 Still one needs to have access to m^-values for which m^L > 3.5. Otherwise the finite- volume effect 
for m n is not under control. 
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Figure 9: Top: combined fit to our nucleon mass (red diamonds) and cr-term data (light green band). 
Shown are the volume-corrected data (determined through the fit) for one volume only. The fit's char- 
acteristics (fit ranges, fixed parameters) are the same as for the fit labeled Sool in Table lRl The inner 
band is the error band for the fit, the outer band illustrates the variation of the fit changing (3 and C3 
by one standard deviations around their phenomenological values [Eq s. l|12| and jl3| ] . For the physical 
point, marked by a (yellow) circle, we use our final estimate for ro [Eq. ( fTop J! Bottom: same as top panel, 
but the combined fit also includes the nucleon mass data of the ETM Collaboration (blue squares) |30| . 
Note that their data is plotted against the 7r + mass. 
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Agreement is also found if one compares with other estimates, for example, with the 
N f = 2 + 1 estimate of the BMW Collaboration [TO], the CSSM results for N f = 2 and 
A// = 2+1 [XTl I2B"1 |3"I] or with one of the ETM estimates for N f = 2 [3D] (see Fig.[8] 
for a comparison). Note that ETM quotes a final result <T p h ys = 64(l)MeV from a fit to 
0(p 3 ) BChPT, however they find also <7 p hy S = 39(12)MeV from a modified fit function 
(see Figs. 13 and 14 of [30 ). Also from our fits to 0(p 3 ) BChPT we could find larger 
values for a, depending on the largest pion mass included in the fits. However, based on 
the systematic procedure we have applied here, we think the value in Eq. ( 16 1 is more 
reliable. 

To compare our results with other recent Nf — 2 calculations we have redone our 
combined fits including the raw (i.e., not corrected for finite-size effects) data of the 
ETM Collaboration (3U] into our analysis. This roughly doubles the amount of data 
which is available at < 300 MeV. Repeating fits for different fixed parameter values 
?3, C2 and C3 we obtain values for tq and <7 p hy S which agree within errors with those of 
Eqs. (15) and (16). The bottom panel of Fig. [9] shows one of these fits: Red diamonds 
represent our data, blue squares ETMC's points, all after subtracting the finite-volume 
corrections. As above, these corrections were determined directly through the fit, using 
points from different volumes but same (r m w ) 2 . For simplicity though, in this figure 
only points from the largest available lattice volumes are shown. The short (green) band 
is from the value rpcr = 0.273(25) |7J which constrains the slope of the fitting function at 



rom v = 0.735. The yellow circle marks the physical point using our rp-value in Eq. (15 



For a comparison, the top panel of Fig. [9] shows the corresponding fit (labeled Sool in 
Table lB~4|) without ETMC's points included. 

We conclude with a note on r^: The value tq = 0.501(10)(11) fm obtained from our 
fits is somewhat surprising. In other studies (including ours) ro was often found to be 
about 0.47 fm (see, e.g., (3UJ l3"2Tf3"4"] ). Here ro was fixed iteratively forcing the fit to be 
self-consistent, i.e. ro = Mjv(ro ■ m p h ys )/Af]J hys . Our value for ro (and that for <7 p hy S ) 
is thus valid as long as the following assumptions are satisfied: (a) O(m^) BChPT is 
sufficient to describe our data below = 435 MeV, (b) discretization effects are really 
negligible and (c) a physical value of 938 MeV is adequate for a two-flavor calculation. 
Whether these assumptions have to be changed or relaxed has to be seen when there 
will be more data available at smaller pion masses and lattice spacings, and also for 
simulations with Nf = 2 + 1 and Nf = 2 + 1 + 1. From the right panel of Fig. [6] however, 
we see a trend for ro: it increases when the upper limit on r^m^ is lowered. 

Note that also a recent study [35] of the kaon decay constant finds a r -value of around 
0.5 fm from a two-flavor lattice calculation. Their estimate is completely consistent with 



our estimate of Eq. ( 15 
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Appendix A. Nucleon mass and pion-nucleon cr-term from BChPT 

In this appendix we derive our fit formulae from a next-to-leading one-loop order (or 
C(p 4 )) calculation of the nucleon mass in covariant baryon chiral perturbation theory 
(BChPT) for two light flavors. The pion-nucleon a term then follows from an applica- 
tion of the Feynman-Hellmann theorem (see Eq.(2)). In addition we need to know the 
connection between the quark mass and the pion mass to 0(p 4 ) [SUSHI. Denoting the 
mass of the degenerate light quarks by me — m u — m& the leading term is given by the 
Gell-Mann-Oakes-Renner relation 

ml = 2B m e . (A.l) 

While the quark mass me and the parameter Bq are scheme and scale dependent, the 
auxiliary variable 

m 2 = 2B m e (A.2) 

is independent of these conventions. 

The generic 0(p ) result for the nucleon mass can be written as 

M N = M + M (1) + M (2) + M (3) + M (4) + 0(p 5 ) . (A.3) 

Evaluating the required loop diagrams with IR-regularization [41] one obtains 

M (1) =0, (A.4) 
M (2) = - 4cim 2 , (A.5) 



uriVi 3(g° A ) 2 m 3 | / m 2 m m m i 

M ' = n, nxo S \ / 1 o arccos 1 log — * } . (A.6) 

167r 2 (F£) 2 ly 4M 2 2M 4M B M 2 ' 



M(i) = ^ ^ + 64^ £ { 8C1 - " - ^ - ^ 2 



3m 4 /c 2 o y A f | (g°) 2 log 



64tt 2 (F°) 2 \ 2 M M 6 M 2 



3ci(^) 2 m 6 I m 2 m arccos 



l 0g - _ _ 2 ^» K . (A.7) 



Here and F® denote the axial coupling constant of the nucleon and the pion decay 
constant in the chiral limit, and Ci, C2, C3 are the standard low-energy constants (see, 
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e.g., [H3). The renormalized counterterm coefficient e£(A) depends on the scale A of 
dimensional regularization in such a way that does not depend on A. 

The generic next-to-leading one-loop result for the sigma term of the nucleon reads 



r (3) 



Utilizing the relation 



a 



dM N (m e ) , dM N (m 2 ) 
mi = to 



dmg 



dm 



one obtains from Eqs.(A.4|-( A. 7) 
<r«=0, 

(2) _ ™2 



r(3) = 



4ciTO 



3(.9^) 2 
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16tt 2 (F0) 
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2M 
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(A.8) 
(A.9) 

(A.10) 
(A.ll) 

(A.12) 



4c 3 



(<£) 2 ' 
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32712(^0)2 



4ci - 2c 3 - 3 
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16tt 2 (F°) 2 M 2 



3 log 



M Q 1Og M 2 



M 2 



3rrT 
4M§ 



rn 



2.U,, 



3/2 



(A.13) 



Finite volume corrections to the nucleon mass can be evaluated from the same Feyn- 
man diagrams. At next-to-leading one-loop order one gets [43j 



ith 



AM< 3 ' 
AM< 4 ' 



AM N (fn 2 ,L) = AM (3) (m 2 ,L) + AM (4) (m 2 ,L) +0(p 5 ) 

3(g0)2M TO 2 
16tt 2 (F0) 2 



dx^2 Kol L ,7 



\n\ sj M 2 x 2 + m 2 (l - x) 



3to 4 x - 
4tt 2 (F0)2 2_, 



AjX^to) fT 2 (£|n|m) 
(2C1 " C3) +C2 (L|n|TO)2 



(A.14) 

(A.15) 
(A.16) 



where Ki is a modified Bessel function. Note that the value n = is omitted in the 
threefold sum over the integers ni, n 2 , n^. 

For the applications in the present paper it is advantageous to consider the nucleon 
mass as a function of the pion mass m^ (the mass of the lowest lying - state in the 
simulation). Therefore we have to convert our expressions for M N {m 2 ) of Eq.([A3| into 
expressions for Mjv^m 2 .). Utilizing the 0(p 4 ) result [55] 
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(A.17) 



we can eliminate the dependence on m at the cost of introducing the renormalized low- 
energy constant Z 3 (A), which cancels the dependence on the renormalization scale A of 
the associated chiral logarithm in Eq. (A. 17). Note, however, that l^(X) can be subsumed 
into an effective coupling eJ(A) via 



e ;(A) = eI(A) + 2Z5(A) 



Cl 



(^°) 2 



(A.18) 



and therefore cannot be determined independently within an analysis of the mass of the 
nucleon. One hnds 
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(A.19) 



Note that we have shifted the couplings g A — > gA, — > i 7 ^ to their physical values, 
consistent to the order at which we are working. This applies also to the finite volume 
corrections, where we may in addition replace m by m ff . 
For the a term as a function of we obtain 
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(A.20) 



Note that in contrast to the chiral extrapolation function of the mass of a nucleon given 
in Eq. (A.19) the dependence on the counter term Z 3 (A) introduced in Eq. ( A.17[ ) cannot 
be absorbed completely into the effective coupling e£(A). In order to calculate a(m,P hvs ) 
we therefore also need information on the numerical size of the coupling Z 3 , which is 
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related to the scale-dependent coupling i?j(A) via 



TO = -g^(/3+log™H • (A.2U 



Finally, expanding the expressions for M^im^) and <t(to w ) given in Eqs. (A. 19 1 and 
(A. 20 1, respectively, in powers of up to 0{rn^) we arrive at the formulae Eqs. Q 



and ([5j) used in our fits. 

Let us now estimate the theoretical uncertainty associated with our next-to-leading 
one-loop BChPT calculation of Mjy due to the truncation of 0(p 5 ) effects (and all higher 
orders). A complete two-loop calculation of the nucleon mass in BChPT, employing a 
reformulated version of infrared regularization, has been worked out in Refs. [23, 44J. 
Truncating the result at O(m^), we find 



M 



N = M + ~k x m\ + k 2 m\ + k 3 m^ log ( -r^ J + fc 4 m^ 

+ k 5 ml log + k 6 ml + 0(m*) . 

The coefficients of this expansion are given by 

k = -4ci, h = ~ 3{9 ° a)2 
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The one-loop expression (coefficients ki - ki) is consistent with Eq.(|4|). In order to 
study the higher order effects we need values of the additional low-energy constants l\ 
and 2d\ 6 - d r w . From [EUDUm] we find l i( X = 0-!38GeV) = 0.027 and (with a consid- 
erable uncertainty) di 6 (A = 0.f38GeV) = — f.76GeV~ 2 . For the scale-independent 
constant d[ s , Ref. derives the value d[ s — — 0.80 GeV~ 2 from the Goldberger- 
Treiman relation. In view of the large uncertainty of d r 1& we use the rough estimate 
2d\§ — di 8 = (—2.5 ± 2.0) GeV~ 2 . This estimate, however, does not reproduce the m w - 
dependence of our nucleon mass data. Using 2d\ e — d\a = 3.0 GeV~ 2 instead results in 
much better fits to the data. The 0(m%) curves in Fig.Fnhave been calculated with both 
these numbers. 

The corresponding expansion of the a term is calculated from the derivative of Mjv 
with respect to the light quark mass according to Eq.Q. Expressed in terms of it 
reads 

a = himl + h 2 ml + h 3 m* log ( — ^ ) + h A ml 

+ h 5 ml log (-^J + h 6 ml + 0{ml) , 
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with the coefficients 



fc - -4c , h- 9(95)2 



h 



3 



64tt(F0)2 ' 

3{g° A ) 2 - 28M ci + 3M c 2 + 12M c 3 
16Af ^ 2 (F°) 2 



~h 4 = 8el-^|-— 3^ ( 3 (s o )2 _ 8MoCi+4MoC3) ^ 



ft 



5 



(i^o) 2 64M tt 2 (^o) 2 
3(^) 2 (40(^) 2 -3) 



10247r 3 (i^; 



4 



S » - 204^(jy { ^ (12a)2 - 1)M ° 

+ 4>r 2 [16 9 » (5(5 - 3i;)Mj + 5(f;) j ( 9 5 - 32M 2 (2< e - <?,))] } . 

Finally, we add a remark concerning the LEC e^. Becher and Leutwyler give an 
estimate for the Delta contribution to ef L (see Appendix D of [45] ) : 

-bl _ 9a ( A 

e l.A " 



6 log UtT" • (A- 24 ) 



48tt 2 A V \ 2M o 

For 5A = 13.0 GeV -1 , A = (1.232 - 0.939) GeV, this amounts to ef^ ~ 7.5 GeV~ 3 . 
Translated to our choice of LECs, this would give 

,- rn , 3(g 2 - 8ciM + c 2 M + 4c 3 Af ) / m phys \ _.„„_ 3 
4e l( A = M„) = 4 61 + 32Mo7r2(f o )2 log J = 7.5 GeV 

=>■ -1 GeV" 3 < e\{\ = 0.138 GeV) < 1 GcV~ 3 , (A.25) 

so this (shifted) LEC may therefore be expected to be of natural size. Such resonance 
saturation estimates usually give only very rough estimates for the higher-order LECs, 
at least in the baryonic sector. Nonetheless, such a range of values for e\ would be 
consistent with the findings of our fits (cf. Appendix B). 



Appendix B. More details on the fits 

Here we give a short summary and discussion of the fit parameters for our (simultane- 
ous) fits to the nucleon mass (and cr-term) data. For both types, fits have been performed 
for different fit ranges and fixed input parameters, c%, c 3 and Z 3 , in order to systematically 
explore the dependence of the fit parameters on this external bias. Table lB~4l summarizes 
the results for the combined fits and Tablc lB~5l those for the stand-alone fits. 

The two tables have to be read as follows: Each line is the result for one particular 
fit. In the first column, a unique key is assigned to each fit. Keys starting with N ("no 



<t") refer to stand-alone fits to the nucleon mass data [i.e., Eq. (10)], while keys starting 



with S ("with er") signify combined fits [Eq. (11)]. The second and third characters are 
either m, o, p or f , depending on the values assigned to the parameters c 3 and i 3 (see 
Scc. l3.3p . If the second (third) character is an o, c 3 (f 3 ) was fixed to c 3 = —4.7 GeV -1 
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Key 


(rom 7r )^ la 


xl 


r M 


ci/ro 


e l/ r 




c zt r o 




r [1/GeV] 


+ Smml 


1.30 


4.470 / 5 


2.26(4) 


-0.34(3) 


-0.142(11) 


1.30 


-2.37 


0.098(17) 


2.53(4) 


* Smol 


1.30 


4.450 / 5 


2.22(4) 


-0.36(3) 


-0.148(10) 


1.32 


-2.39 


0.104(17) 


2.51(4) 


* Smpl 


1.30 


4.810 / 5 


2.19(4) 


-0.39(3) 


-0.155(10) 


1.33 


-2.42 


0.111(18) 


2.48(4) 


* Soml 


1.30 


4.550 / 5 


2.31(4) 


-0.28(3) 


-0.022(11) 


1.28 


-1.83 


0.084(18) 


2.57(4) 


* Sool 


1.30 


4.060 / 5 


2.29(4) 


-0.29(3) 


-0.024(11) 


1.29 


-1.84 


0.089(17) 


2.56(4) 


* Sopl 


1.30 


3.690 / 5 


2.26(4) 


-0.31(3) 


-0.027(10) 


1.30 


-1.85 


0.094(17) 


2.53(4) 


Spml 


1.30 


6.990 / 5 


2.37(4) 


-0.22(3) 


0.087(12) 


1.26 


-1.30 


0.070(22) 


2.61(5) 


* Spol 


1.30 


6.380 / 5 


2.35(4) 


-0.23(3) 


0.087(11) 


1.27 


-1.31 


0.074(22) 


2.60(5) 


* Sppl 


1.30 


5.770 / 5 


2.33(4) 


-0.25(3) 


0.087(11) 


1.28 


-1.32 


0.079(22) 


2.58(5) 




1.30 




2.28(4) 


-0.31(3) 


-0.048(11) 






093f 1 H) 


2.54(4) 


* Sf ml 


1.30 


4.180 / 4 


2.28(7) 


-0.31(7) 


-0.089(127) 


1.29 


-2.13(58) 


0.092(30) 


2.55(7) 


* Sfol 


1.30 


3.850 / 4 


2.26(7) 


-0.32(7) 


-0.073(125) 


1.30 


-2.06(56) 


0.095(30) 


2.54(7) 


* Sfpl 


1.30 


3.620 / 4 


2.25(8) 


-0.33(7) 


-0.054(122) 


1.31 


-1.98(55) 


0.098(30) 


2.52(7) 


avc. 


1.30 


Xr < 1-3 


2.26(7) 


-0.32(7) 


-0.072(125) 


— 


-2.05(56) 


0.095(3) 


2.54(7) 


* Smm2 


1.60 


5.140 / 7 


2.25(3) 


-0.35(2) 


-0.146(5) 


1.30 


-2.37 


0.100(9) 


2.53(3) 


* Smo2 


1.60 


5.100 / 7 


2.22(3) 


-0.36(2) 


-0.150(5) 


1.32 


-2.39 


0.105(9) 


2.51(3) 


* Smp2 


1.60 


5.460 / 7 


2.19(4) 


-0.38(2) 


-0.154(5) 


1.33 


-2.42 


0.110(10) 


2.48(3) 


* Som2 


1.60 


5.270 / 7 


2.30(3) 


-0.29(2) 


-0.028(5) 


1.29 


-1.83 


0.089(10) 


2.56(3) 


* Soo2 


1.60 


4.730 / 7 


2.28(3) 


-0.31(2) 


-0.029(5) 


1.30 


-1.84 


0.093(9) 


2.55(3) 


* Sop2 


1.60 


4.350 / 7 


2.25(4) 


-0.32(2) 


-0.030(5) 


1.30 


-1.86 


0.097(9) 


2.53(3) 


* Spm2 


1.60 


7.820 / 7 


2.35(3) 


-0.24(2) 


0.079(5) 


1.27 


-1.31 


0.078(12) 


2.60(4) 


* Spo2 


1.60 


7.100 / 7 


2.33(3) 


-0.25(2) 


0.080(5) 


1.28 


-1.32 


0.081(12) 


2.59(4) 


* Spp2 


1.60 


6.430 / 7 


2.31(4) 


-0.26(2) 


0.081(5) 


1.28 


-1.32 


0.084(11) 


2.57(3) 


avc. 


1.60 


Xr < 1-3 


2.28(4) 


-0.31(2) 


-0.033(5) 






0.095(10) 


2.54(3) 


* Sfm2 


1.60 


4.870 / 6 


2.27(6) 


-0.32(6) 


-0.095(123) 


1.30 


-2.15(57) 


0.096(18) 


2.54(5) 


* Sfo2 


1.60 


4.510 / 6 


2.26(6) 


-0.33(6) 


-0.077(121) 


1.30 


-2.07(55) 


0.098(18) 


2.53(5) 


* Sfp2 


1.60 


4.280 / 6 


2.24(6) 


-0.33(6) 


-0.055(118) 


1.31 


-1.97(54) 


0.099(18) 


2.52(5) 


avc. 


1.60 


Xr < 1-3 


2.26(6) 


-0.33(6) 


-0.08(12) 




-2.1(6) 


0.098(18) 


2.53(5) 


Smm3 


3.00 


50.650 / 9 


2.11(3) 


-0.48(1) 


-0.178(5) 


1.36 


-2.47 


0.144(11) 


2.43(7) 


Smo3 


3.00 


44.340 / 9 


2.08(3) 


-0.49(1) 


-0.180(5) 


1.38 


-2.50 


0.145(10) 


2.40(7) 


Smp3 


3.00 


38.890 / 9 


2.05(3) 


-0.51(1) 


-0.182(6) 


1.39 


-2.53 


0.146(9) 


2.37(6) 


Som3 


3.00 


36.360 / 9 


2.18(3) 


-0.40(1) 


-0.045(5) 


1.33 


-1.89 


0.127(10) 


2.49(6) 


Soo3 


3.00 


31.150 / 9 


2.16(3) 


-0.40(1) 


-0.043(5) 


1.34 


-1.91 


0.128(9) 


2.47(6) 


Sop3 


3.00 


26.420 / 9 


2.14(3) 


-0.41(1) 


-0.041(5) 


1.35 


-1.92 


0.129(8) 


2.45(5) 


Spm3 


3.00 


29.170 / 9 


2.25(3) 


-0.32(1) 


0.072(4) 


1.30 


-1.34 


0.110(9) 


2.54(5) 


Spo3 


3.00 


25.220 / 9 


2.24(3) 


-0.33(1) 


0.075(4) 


1.31 


-1.35 


0.111(9) 


2.52(5) 


Spp3 


3.00 


21.550 / 9 


2.23(3) 


-0.34(1) 


0.078(5) 


1.31 


-1.35 


0.112(8) 


2.51(5) 


Sfm3 


3.00 


27.120 / 8 


2.33(9) 


-0.24(9) 


0.203(159) 


1.27 


-0.68(75) 


0.090(61) 


2.59(13) 


Sfo3 


3.00 


24.190 / 8 


2.31(10) 


-0.26(10) 


0.180(175) 


1.28 


-0.82(82) 


0.095(64) 


2.57(14) 


Sfp3 


3.00 


21.180 / 8 


2.28(12) 


-0.28(11) 


0.159(188) 


1.29 


-0.95(87) 


0.099(65) 


2.55(15) 



Tabic B.4: Parameters from our combined fits to the nucleon mass and cr-term data. Keys in the first 
column indicate the features of the fits. All keys start with S . . indicating that the cr-term data has been 
included when fitting. Entries are ordered according to increasing (roWilLxi as specified in column 2. 
The x 2 -value together with the number of degrees of freedom (ndf) is given in column 3. Columns 4 to 
8 list the fit parameters Mo, ci, eJ(0.138MeV), C2 and C3 in units of r<j. The corresponding estimates 
for rrj(7 at the physical point and for rrj are given in column 9 and 10. Since r<j is not a fit parameter, 
but has been iteratively fixed as explained in the text, the error for r<j is that of M^(tq ■ m^ hyB ) / M^ yB . 
Lines tagged with a * indicate fits where x 2 / n 4f < 1-3. Parameters of these fits enter the weighted 
averages given in the lines starting with "ave. ". 



22 




Figure B.10: Volume dependence of our nucleon mass data at (/3,k) = (5.29,0.13632) (left) and 
(5.40,0.13660) (right) resulting from different combined fits, labeled Smpl, Sool and Spml in Table lB~4l 
They illustrate the maximum variation of the finite-volume corrections for all our combined fits [for 
(roiriTr) 2 < 1.3] with the parameters C3 and £3. 



Key 


( r m «)L 


x xl 


r M 


ci/r 




C2/r 


C3A0 


ro Cphys 


r [1/GeV] 


Nmol 
Nool 
Npol 


1.30 
1.30 
1.30 


4.760 / 4 
3.510 / 4 
4.450 / 4 


2.19(9) 
2.23(9) 
2.26(9) 


-0.39(6) 
-0.34(6) 
-0.29(6) 


-0.155(11) 
-0.030(11) 
0.088(11) 


1.33 
1.32 
1.31 


-2.42 
-1.88 
-1.35 


0.111(40) 
0.103(35) 
0.096(40) 


2.48(10) 

2.50(9) 

2.53(10) 


avc. 


1.30 


%l < 1-3 


2.23(9) 


-0.34(6) 


-0.032(11) 






0.10(4) 


2.50(10) 


Nlol 


1.30 


3.510 / 3 


2.23(10) 


-0.34(8) 


-0.024(133) 


1.32 


-1.85(58) 


0.103(50) 


2.51(12) 


Nmo2 
Noo2 
Npo2 


1.60 
1.60 
1.60 


5.280 / 6 
4.200 / 6 
5.170 / 6 


2.20(7) 
2.23(7) 
2.26(7) 


-0.38(4) 
-0.34(4) 
-0.30(4) 


-0.153(8) 
-0.031(8) 
0.086(7) 


1.32 
1.32 
1.31 


-2.41 
-1.87 
-1.35 


0.109(23) 
0.103(21) 
0.097(24) 


2.49(7) 
2.51(6) 
2.52(7) 


avc. 


1.60 


xl < 1-3 


2.23(7) 


-0.34(4) 


-0.022(8) 






0.103(22) 


2.51(6) 


Nfo2 


1.60 


4.200 / 5 


2.23(8) 


-0.34(6) 


-0.029(132) 


1.32 


-1.87(58) 


0.103(26) 


2.51(7) 


Nmo3 
Noo3 
Npo3 
Nfo3 


3.00 
3.00 
3.00 
3.00 


35.320 / 8 
17.010 / 8 
11.140 / 8 
10.810 / 7 


1.92(4) 
2.04(4) 
2.14(4) 
2.18(11) 


-0.57(1) 
-0.46(1) 
-0.37(1) 
-0.34(10) 


-0.191(9) 
-0.030(8) 
0.099(7) 
0.150(156) 


1.47 
1.41 
1.36 
1.34 


-2.68 
-2.00 
-1.40 
-1.15(69) 


0.156(11) 
0.138(9) 
0.121(8) 
0.114(43) 


2.24(8) 
2.35(6) 
2.43(4) 
2.46(12) 



Table B.5: Parameters from our fits to the nucleon mass alone. All keys start with W. . indicating 
that no CT-term data has been included. The remaining character sequence has the same meaning as in 
TablclB~4l 
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(Is = 3.2); if instead this character is p (m), C3 and I3 were fixed to these values, plus 
(minus) one standard deviation. If the second character is f , C3 was not fixed but left as 
a free fit parameter. Since C2 is known to a much better precision than C3 [see Eq. ( 13 1] 
it is not varied, but always set to its phenomenological value C2 = 3.3 GeV -1 . Each key 
ends with an integer labeling, in increasing order, the upper limit ('"oW 7r )j lax of the fit 
interval 

where (rom T )^ ax is either 1.3, 1.6 or 3.0 (given in the second column). In all cases we 
require L/ro > 3. To reduce the finite- volume effect in the pion mass, we exclude data 
points for all those (/?, k) combinations for which there is not at least one value for rom,,- 
available that satisfies m^L > 3.5. 



All (fixed and floating) fit parameters listed in Tables B.4 and B.5 are given in units 



of tq. The corresponding self-consistent physical value of tq reached itcrativcly for each 
fit (see Eq. ([£])) — is given in the last column in GeV -1 . That is, if one is interested in the 
physical value for a particular fit parameter, this parameter has to be multiplied by the 
corresponding power of Tq given at the end of the same line. In TableE! we also give 
weighted averages of the fit parameters for fits with a x 2 - < 1-3 (see lines starting with 
"ave."). 

Let us now comment our fits. The fact that our fits with (rom^) 2 < 1.6 yield (in 
the majority of cases) Xr _vames around one indicates that the constraint L/ro > 3.0 on 
the spatial lattice extension has been sufficient to correct for the finite-volume effect in 
the data. This can be seen also from Fig. |B.10"} where we compare data for r M^ at 
fixed (/?, k) but different L/r to the fitted function at the corresponding r m n . If we 
had relaxed the constraint L/ro > 3, the finite- volume effect could not be completely 
compensated for data points with L/ro < 3. This is, for example, the case for our points 
at (P,k) = (5.40,0.13640) from a 24 3 x 48 lattice. 

Looking at the fit parameter e£, our fits are less robust, however. In fact, the sign 
of e\ is strongly correlated with the value we choose for C3. If we set, for example, 
C3 = —4.7 + 1.3 GeV -1 all fits result in a positive e\, while, if we set C3 = —4.7 GeV -1 or 
C3 = —4.7—1.3 GeV -1 , e\ comes always out negative. The value we choose for I3 does not 
affect the sign of ej, yet I3 has a minor effect on e"s absolute value. To our knowledge, 
nothing is really known about the sign of e\, so we do not restrict it. From our data we 
can say it tends to be negative if C3 < — 4.7GeV and positive if c 3 > —3.4 GeV -1 . If we 
leave C3 as free fit parameter, we obtain values around —5.0 GeV -1 for C3 (±1.5GeV _1 
statistical uncertainty), and e 1 is consistent with zero within errors. This is found for 
fits where (ro"i,)J,„ is either 1.6 or 1.3. Interestingly, the range of fitted values for e\ 
lies in the ballpark expected from BChPT (see Eq. (A.25l). 

A correlation with c 3 is also seen for ro<7 p hy S and r (see columns 9 and 10 in Tables 
and B.5). For larger, i.e., less negative values of C3, the results for ro<x p hys tend to 



B.4 



smaller numbers, while ro tends to larger values. 

For the reader's convenience, we visualize the variation of the values for ro and rotT p hy S 
in Fig. |B.ll] (only for results from Table lB~4|) . From top to bottom, panels are ordered 
with decreasing (^o m 7r)maxi while within each panel, points are grouped according to the 
values chosen for C3 and Z3. Symbols distinguish different c 3 , neighboring points with the 
same symbol are for different Z3. The color intensity of each point is related to the Xr~ 
value of each point. When applicable each panel also shows the corresponding weighted 
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average. 
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